Decoherence and relaxation in the interacting quantum dot system 
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In this paper we study the low temperature kinetics ol the electrons in the system composed of a 
quantum dot connected to two leads by solving the equation of motion. The decoherence and the 
relaxation of the system caused by the gate voltage noise and the electron-phonon scattering are 
investigated. In order to take account of the strong correlation of the electrons in this system, the 
quasi-exact wave functions are calculated using an improved matrix product states algorithm. This 
algorithm enables us to calculate the wave functions of the ground state and the low lying excited 
states with satisfied accuracy and thus enables us to study the kinetics of the system more efi'ectively. 
It is found that although both of these two mechanisms are proportional to the electron number 
operator in the dot, the kinetics are quite different. The noise induced decoherence is much more 
effective than the energy relaxation, while the energy relaxation and decoherence time are of the 
same order for the electron-phonon scattering. Moreover, the noise induced decoherence increases 
with the lowering of the dot level, but the relaxation and decoherence due to the electron-phonon 
scattering decrease. 

PACS numbers; 71.27.-|-a,72.15.Qm,73.21.La 
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I. INTRODUCTION 

In the past two decades there have been many ex- 
perimental and theoretical investigations on the artifi- 
cial strong interacting quantum system, especially af- 
ter the observation of Kondo effect in the quasi one- 
dimensional (ID) system composed by a semiconductor 
quantum dot (QD) connected to two leads. The Kondo 
effect of this system produces a perfect transparency 
with unitary conductance in the symmetric configura- 
tion £o — ~U/2, with £q and U being the QD level and 
on-site Coulomb repulsionjii^i^ii^ Most of the studies fo- 
cus on the static transport properties. In recent years, 
the kinetics and time evolution of these systems also at- 
tracted much attentionj^iiiSii^ as the study of the kinetics 
and time evolution enables one to attack the strong corre- 
lated system from different perspectives. Different meth- 
ods have been employed to study the response of the sys- 
tems to the ac- modulated or step switching bias and/or 
gate volt age i^iiSiiiii^ Theoretically, the temporal evolu- 
tion of Kondo-like systems has been studied by using 
noncrossing approximation )^i^°'^^ quantum Monte Carlo 
simulatiouflS time-dependent numerical renormalization 
groupfi^ and time-dependent density-matrix renormal- 
ization group (DMRG) methods. However, only few 
studies have considered real dissipation mechanisms, 
such as noise and the electron-acoustic (AC) phonon in- 
teraction which dominate the decoherence and relaxation 
at low temperature^^ 

In this paper, we study the time evolution of the QD 
system under the influence of the gate noise and acous- 
tic phonon by solving the equation of motion. Since the 
electrons in this kind of system are mutually correlated, 
an improved version of the matrix product states (MPS) 
algorith m "'^^i^^'^^i^^'^^i^^ is employed to study the electron 
system in order to keep the many-body effect of the sys- 
tem accurate. The paper is organized as following. In 



the second section, we present the system Hamiltonian 
and the equation of motion for the time evolution of the 
system in the presence of the gate noise and electron 
phonon interaction. In the third section we give a brief 
introduction to the MPS algorithm and discuss how it 
can be further improved to study the kinetics of the QD 
system. We show the numerical results of the decoher- 
ence and relaxation in the fourth section and summarize 
in the last section. 



II. SYSTEM MODEL 

The system we study is composed of a semiconductor 
QD connected to two leads. The Hamiltonian is written 



H = H, + H, 



ph 



(1) 



Here He and Hph are the Hamiltonian of electron and 
phonon, Hgp describes the electron-phonon interaction. 
The electron Hamiltonian of a iV-site QD system is 



Uno^nol + UCla+ia + h.c. , (2) 



with Cjg.{Cia) being the creation (annihilation) opera- 
tor of electron with spin cr (=ti) at the site i. The sites 
with index i = 1, 2, • • • , [N/2]{-[N/2], • • • , -2, -1) are 
the sites on right (left) lead, while site with i = is for 
the QD. In this paper, the on-site energy of leads is set 
to be zero and = t for z = • • • , —3, —2, 1, 2, • • • . U de- 
scribes the Coulomb repulsion in the QD. The QD level Eq 
can be controlled by gate voltage and is subjected to the 
influence of the noise. For the phonon related part, we 
only consider the electron interaction with background 
longitudinal AC phonon via deformation potential. The 
phonon Hamiltonian and the electron-phonon scattering 
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III. NUMERICAL SCHEME 



and 



qA 



Hep^Y^Y. ^OaM^ibq + bl^) , (4) 



respectively. &q(f'q) is the creation (annihilation) opera- 
tor of longitudinal AC phonon with wave- vector q, whose 
frequency is Wq = Vsiq with Vsi standing for the lon- 
gitudinal sound velocity. Mq (x ^(^'o|e''i'"|5'o) is the 
corresponding coupling matrix. Here is the electron 
wave- function in the dot. It is chosen to have Gaussian 
form exp(— r^/2a^) for simplicity, with a denoting the 
diameter of the QD. 

The kinetics of the system is studied by the temporal 
evolution of the density matrix /o(t), whose diagonal ele- 
ments pn,n{T) and off-diagonal ones pn,m{T) stand for the 
population of \n) state and the coherence between \ n) and 
\m) states at time r respectively. Here \n) and |m) are 
the eigenstates of electron Hamiltonian free from the gate 
noise and electron-phonon interaction. Using Markovian 
approximation and secuUar approximation, one can write 
down the equation of the motion for the density matrixi^^ 

^ = ~i[H,,p] - ^[i^Ls, p[t)] J2 

€ 

nl{£)piT) + p{T)no{e)nl,{e) ~2nl{e)p{T)no{e)] . (5) 

Here Hls ~ A(e)no(e)no(e) is the energy shifting 
due to electron-phonon interaction. A(£) and r(e) are 
the real and imaginary parts of the following formula 

/•OO 

A{e) + trie) = l^-^ql' / dTe''^D{(i, r) , (6) 

q "'O 

with Z)(q, r) being the phonon Green function. no(e) ~ 
J2n m^s,En-E,„\n'){n\no\ni){m\. For the strong electron- 
optical phonon interaction, Hj^g plays important role in 
the system properties;^ However, for the weak electron- 
AG phonon interaction we consider in this paper, it only 
slightly rcnormalizes the QD level and the Gounlomb in- 
teraction, which is hard to detect experimentally. There- 
fore this shifting is simply omitted in this paper. For 
electron-AG phonon interaction via deformation poten- 
tial r(e) takes the following form 

r(e) = Ap|en(l + nB{e))9{e) + nB{-e)e{-e)] , (7) 



where n_B(e) = !/{• 



,e/T _ 



1) is the Bose function at tem- 



perature T, and 9{e) is the Heaviside step function. Ap 
depends on the material and structure dependent param- 
eter as well as on the electron-phonon coupling. 



To solve the equation of motion, one needs the the 
wave- function of eigenstates of the electron system to ob- 
tain the matrix elements of density matrix and hq. Since 
the electrons governed by the Hamiltonian [Eq. ^] form 
a mutually correlatted system, it is important to take the 
correlation into account in the study of the kinetics at low 
temperature. Here we use an improved MPS algorithm to 
obtain the many-body wave-functions so that the strong 
correlation of the electron system can be automatically 
taken into account. 

For ID ground state problem, MPS is known to be 
equivalent to BMRG^^^^^ DMRG was first proposed to 
study the ground and low-lying excited states of quan- 
tum systems and later was further extended to 
the simulation of the time evolution, calculation of exci- 
tation spectra and finite temperature properties of quan- 



tum systems 
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The success of DMRG is eventually 



understood by its connection to MPS J^i^^i^ In the MPS 

algorithm, the wave-function of a iV-site lattice is rep- 
resented by a group of A-matrices whose dimension is 
usually much smaller than the dimension of Hilbert space 



J2 Tt{A'''---A'^-'A'^A''^+'---A'''} 



\Si - ■■ Si-iSiSi+i ■ ■ ■ sl) 



(8) 



with Si representing the local state index at i-th site. 
To start the calculation, an initial wave- function is given 
by other calculations or randomly generated. The wave- 
function is then gradually optimized by a process called 
"sweep" : The wave-function is optimized by only mini- 
mizing the energy with regard to the local A-matrix at 
one or two "center" sites while keeping other A-matrices 
unchanged. One then moves the "center" site to the left 
or right neighbor by targeting to the ground state2ii2^ 
or performing singular value decomposition.^^ The pro- 
cess repeats until a desirable accuracy is achieved or the 
accuracy can not be further improved by more sweeps. 
In traditional DMRG algorithm, if more than one eigen- 
state are to be calculated, they should be targeted simul- 
taneously during the sweep process. This requires larger 
A-matrix dimension and more computing resources. The 
more eigen-states are required, the less accurate they are 
for a fixed A-matrix dimension. Therefore, the tradi- 
tional DMRG algorithm is usually limited to the calcula- 
tion of the ground state and few low lying excited states. 
Here we propose an algorithm to improve the calculation 
of the excited states by using MPS. 

In our algorithm, instead of calculating all of the re- 
quired eigen-states simultaneously, they are calculated 
step by step. We first calculate the ground state by us- 
ing the traditional DMRG method and obtain the wave- 
function in the normal MPS form like Eq. dS]). To ob- 
tain the first excited state, we again generate an initial 
wave-function then gradually optimize it through sweep. 
However, during the sweep, we first orthogonalize the ex- 
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FIG. 1: Relative error of the ground state energy vs. sweep 
step for the system with N = 17, and Ni; ^ = 8. The 
hopping between the leads and the QD is O.lt and the on-site 
interaction U = 6t. The filled dots are the results of improved 
algorithm; The open boxes are those of single sweep with 
randomly generated initial wave-function. The curves are the 
guilds for the eyes. 

cited state and the ground state by performing Gramni- 
Schmidt orthogonalization on the local A-matrix at the 
center site then minimize the energy for the excited state. 
In this way, we can obtain the wave-function of the first 
excited state in the normal MPS form when the sweep 
converges. Since the wave-function of ground state is in 
the normal MPS form, the sweep of the first excited state 
does not change the ground statci^ Similarly, to calcu- 
late the higher excited states, we first perform the or- 
thogonalization of wave-function to all of the lower states 
then optimize the energy during the sweep. By repeating 
this process, we can get the excited states step by step 
without sacrificing the accuracy of the lower states. 

In our algorithm, the normal MPS form of wave- 
functions are required to carry out the sweep of the ex- 
cited states. There are two major algorithms for sweep: 
two-site sweep and one-site sweep i^^i^^ Two-site sweep 
gives better eigen-energy since it uses a larger Hilbert 
space and therefore it is more memory hungry and CPU 
time consuming. Moreover, the wave-function obtained 
by two-site sweep is no longer in the normal MPS form. 
For example, if the center sites are i — 1 and i, the cor- 
responding matrix of the wave-function are ••• 
j^si-i.si^ • • • . One has to break the matrix A^'^i'*' 

into the multiple of two Matrices with larger dimension 
to get the normal MPS form. The wave-function loses 
its optimality if these two matrix are truncated to the 
original dimension. On the other hand, one-site sweep 
is cheaper to carry out. More importantly, the wave- 
function of one-site sweep algorithm is always in MPS 
form, therefore it keeps the optimality when the center 
site moves. However, one-site sweep can easily fall into 



some local optimal points and usually fails to give satis- 
fied eigen-energy and wave-function unless a good initial 
wave- function is given to start the sweep. Our solution 
to the dilemma is to use the combination of two-site and 
one-site sweep. That is, we first use two-site sweep to im- 
prove a randomly generated wave-function and truncate 
the improved wave-function to the normal MPS form and 
use it as the initial wave-function for the one-site sweep 
algorithm. It is expect that, the truncated wave-function 
would be a much better initial function even when the 
initial two-site sweep is not converged. 

In the following we first use this hybrid sweep algo- 
rithm to study the ground state of the QD system to 
demonstrate its feasibility. In Fig. [1] the relative error of 
ground state energy is plotted as a function of the num- 
ber of the sweep. In the calculation, we keep the most 
relevant 256 states, and the "true" ground state energy 
Eg is obtained by two-site sweep with 400 states kept. 
In the calculation, the total site number, the spin- up and 
-down electron number are N = 17, and N-^ ^ Ni = 8 
respectively. The hopping between the QD and the leads 
is chosen to be t-i = to = Q.lt. In the QD, the on-site 
energy Eq = —2t and the Coulomb repulsion U = 6t. 
For comparison, we also plot the result of pure one-site 
sweep with a random initial wave function. It is noted 
that the first three steps of our new algorithm are two- 
site sweep. One can see from the figure that the initial 
two-site sweep gives very good start point for one-site 
sweep even when the initial two-site sweep is not con- 
verged. The ground state energy obtained by this al- 
gorithm is very close to that obtained by the two-site 
sweep. More importantly, with this hybrid sweep algo- 
rithm, one not only gets an accurate eigen-energy but 
also a very accurate wave-function with cheaper price. 
This sweep algorithm also works for the excited states. 
For higher excited states, the error in the eigen-energy 
is usually larger than those of ground state and lower 
excited states. However, when we keep up to 256 states 
during the sweep the accuracies of the lowest 10 states 
are desirable. 

We now use this method to calculate the wave- 
functions and the eigen-energies of the system at different 
gate voltages. The ground state energy and the elec- 
tron occupation number in the QD of the system with 
N ^ 9,Nf = Ni = 4: are plotted in the Fig. |1 One can 
see that for positive on-site energy, the electron occupa- 
tion in the QD is small and the ground state energy is 
almost independent of gate voltage. For large negative 
QD level, the QD is doubly occupied and the ground 
state energy become a linear function of the QD level. 
The slope of the function is 2, which agrees with the re- 
sult of iVo ~ 2. In the regime of —6t < Eq < —2t, roughly 
corresponding to the transparent regime {—U<eo<0), 
it can be seen that the rate of the changing of occupation 
vs. the changing of the QD level is distinctly slower than 
that of its neighbor. The abnormality in this regime may 
be related to formation of Kondo singlet which is the su- 
perposition of the localized and the delocalized states. 
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FIG. 2: Ground state energy Eg (solid curve) and the occu- 
pation A^o (dashed curve) of the QD vs. the QD level. Noted 
that the axis for A'o is on the right. The hopping between the 
leads and the QD is O.lt and the on-site interaction U = 6t. 



The site and the electron numbers are N = 9, N-^ 
respectively. 



Since the delocalized states have lower occupation num- 
ber in the QD, the formation of the Kondo singlet slows 
down the changing rate of the occupation. 



IV. DEPHASING AND DECOHERENCE 

Once we have the wave- functions, we can simulate the 
kinetics of the interacting system under the gate volt- 
age noise and the electron-phonon interaction. In the 
low temperature regime, it is expected that only the low- 
est few states are involved. Therefore we can study the 
kinetics by solving the equation of motion in the trun- 
cated Hilbert space composed of lowest ten states. We 
first study the decoherence caused by the noise on the 
QD level by solving the equation of motion. In the pres- 
ence of the noise, the QD level becomes Eq = Eq + V{t), 
where Eq is the average QD level, and V{t) is the noise 
due to the fluctuation of gate voltage. For a system 
free from noise, the coherence between different states, 
i.e. the off-diagonal elements of density matrix, evolve as 
/5n.m(0)e~*'-^"~^™-''^ whose amplitudes do not decay with 
time. The noise irregularly shifts the energy levels with 
different amounts for different states and causes transi- 
tion between the states. Since the amplitude of the noise 
randomly changes, different states eventually lose their 
phase information over time. As a result the amplitude of 
the coherence decays with time. In the Fig. [3] (a) we plot 
the time evolution of the coherence between the ground 
state and all of the excited states J2e \Pge{''')\ under dif- 
ferent conditions. To simulate the kinetics of suddenly 
switch of gate voltage, the initial wave-function wave- 
function is chosen to be the unperturbed wave function of 
zero gate voltage in the truncated space. The density ma- 
trix is averaged over 1000 samplings with the noise V{t) 
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FIG. 3: a) The temporal evolution of the averaged coherence 
between the ground state and the excited states for the QD 
level eo = —3t, AV — 0.5t and r„ — 40 /t. b) The decoherence 
time as a function of noise characteristic time r„ for different 
situations: Solid curves are those for AV = 0.5t and dashed 
ones for AV = It; Filled circles and open circles are for eo = 
—3 and —2t respectively. 



evenly distributed over the regime [— AV, AV] to reflect 
the randomness of the noise. One can see from the figure 
that the coherence damped oscillates with time. The en- 
velope of the temporal evolution curve can be fitted by an 
exponential function cx e~'^^'^'' with a decoherence time 
Trf. The decoherence times under different conditions are 
plotted in the Fig. [3] (b). One can see that decoherence 
time decreases with the increase of the noise amplitude 
and the noise characteristic time t„. The on-site energy 
of the QD also affect the decoherence rate, the higher QD 
level the slower the system loses its coherence. The boost 
of the decoherence rate by the increase of noise level is 
quite easy to understand. When the amplitude of the 
noise increases, the transition rates and the difference of 
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energy shifts become larger. Therefore different states 
lose the phase information more quickly. The decrease of 
the decoherence rate with the decrease of t„ is due to the 
suppression of the unitary time evolution by the irregular 
change of the noise. This is similar to the motion narrow- 
ing effect)^ From the Hamiltonian, one can see that the 
perturbation of the noise is proportional to the electron 
number in the QD. Therefore, the average QD level also 
affects the kinetics of the QD system. The lower the QD 
level, the larger the electron number in the QD, and the 
smaller the decoherence time. It should be noted that 
although the gate noise is efficient to remove the phase 
coherence, it is very inefficient to bring the system to the 
thermal equilibrium state as it does not directly carry the 
energy away from the electron system. The relaxation to 
the thermal equilibrium caused by the noise is more than 
ten times slower than the decoherence. 



laxation and decoherence decrease with the lowing of the 
QD level. The differences rise from the fact that electron- 
phonon scattering is inelastic scattering. The electron 
system approaches to thermal equilibrium through the 
absorbing or emitting of phonon. As the electron sys- 
tem loses both energy and phase information by absorb- 
ing/emitting phonon, r,. and are of the same order. 
Furthermore, the main contribution of no operator to 
the phonon induced relaxation/decoherence comes from 
the off-diagonal elements no(£) with £ ^ 0. For lower 
QD level, the electrons are more deeply trapped in the 
QD and it is harder for them to hop to higher-energy 
states. On the other hand, the main contribution of no 
to the noise induced decoherence is the diagonal terms, 
that is the average electron occupation number in the 
QD. Therefore the lower QD level, the quicker the noise 
induced decoherence. 



120 - 



40 - 



G--0 




J I I I I I I I I l_ 



0.2 



0.4 0.6 

Tit) 



0.8 



FIG. 4: The relaxation time (solid curves) and the deco- 
herence time (dashed curves) as functions of temperature T. 
Filled and open circles are for eo = —3 and —2t respectively. 

We now study the relaxation and decoherence due to 
the electron-phonon interaction. In the Fig. 3] we plot 
the relaxation time Tr and the decoherence time Td as 
functions of temperature T under different gate volt- 
ages. Noted that Ap is chosen to be t in the calcula- 
tion. One can see from the figure that both and 
increases as the temperature decreases and saturate at 
low temperature as the electron-phonon scattering sat- 
urates to the emitting phonon limit F(e) ~ Ape^9{e). 
There are two qualitative differences between phonon 
induced relaxation/decoherence and noise induced re- 
laxation/decoherence, even through their perturbation 
Hamiltonians arc both proportional to no, the electron 
number operator in the QD. Firstly, the phonon induced 
relaxation is a little bit faster than the decoherence, while 
noise induced relaxation is more than ten times slower 
than the decoherence. Moreover, the phonon induced re- 



V. CONCLUSION 

In conclusion, we propose an improved matrix product 
states algorithm to calculate the excited states and use 
this algorithm to study the relaxation and decoherence 
caused by the noise and electron-phonon interaction in 
the interacting QD system. Although both of these two 
mechanisms are proportional to the electron number op- 
erator no in the QD, the kinetics due to these two mecha- 
nisms are quite different. The noise shifts the energy with 
different amounts for different states and causes transi- 
tion between states. The irregular change of the noise 
results in the losing of the phase information of the elec- 
tron system and hence of the decoherence. However, the 
noise does not directly carry energy from the electron 
system, therefore the energy relaxation due to noise is 
very inefficient. When the electron-phonon interaction is 
present, the electron system relaxes to the thermal equi- 
librium by absorbing or emitting phonon. As the elec- 
tron system loses both energy and phase information af- 
ter the scattering, the energy relaxation and decoherence 
time are the same order for electron-phonon scattering. 
Moreover, the main effect of the noise comes from the 
diagonal terms of no, while for electron-phonon scatter- 
ing it is from the off-diagonal terms. As a result, the 
noise induced decoherence increases with the lowering of 
the QD level, but the relaxation and decoherence due to 
electron-phonon scattering decrease. 
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